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Abstract 

We give a systematic way to compute higher orders in the l/mh expansion in inclusive 
semi-leptonic decays at tree level. We reproduce the known 1/m^ terms and compute the 
1 /mf terms at tree level. The appearing non-perturbative parameters and the impact on 
the determination of Vet are discussed. 



1 Introduction 



Inclusive semi-leptonic decays are the cleanest way to access the matrix elements of the CKM 
matrix. In order to achieve precision in the determination of the CKM parameters a reliable 
theoretical framework is needed in addition to precise data. Prom the off-diagonal CKM matrix 
elements only V^s and Vet arc known very precisely from direct measurements [1], and both 
determinations can be related to a clean theoretical treatment: while the theoretical machinery 
for the determination of Vus from K iriue decays is chiral perturbation theory, the theoretical 
basis for the precision determination of Vcb is the Heavy Quark Expansion (HQE). 

The HQE [2, 3, 4, 5] has become a very reliable theoretical tool, in particular for semi- 
Icptonic decays. At present, data on inclusive semi-leptonic decays B X£i?£ are so precise 
that not only a determination of Vcb at the level of 2% accuracy is possible [6, 7, 8, 9], but 
also a check of the consistency of the HQE can be performed by determining the parameters 
of the HQE in different ways [10] . These parameters are the kinetic energy parameter /i^^ and 
the chromo- magnetic moment /iq at order 1/m^ and the Darwin term po and the spin-orbit 
term pLS at order We shall stick here with the kinetic scheme, the alternative IS" scheme 

[8, 9] yields similar precision. Using the moments of the hadronic invariant mass spectrum and 
the charged lepton energy spectrum these parameters can be consistently determined with an 
accuracy of about ten percent. 

Within the HQE the order of the moments is related to the order in the l/m?, expansion 
[11, 12]. For example, the moments ((mf, — 2E^)'^) for the case of charmless semi-leptonic 
B decays are determined by the contributions of the order Thus in order to exploit 

precise measurements of the lepton energy spectrum, i.e. measurements of higher moments, 
it is mandatory to perform the theoretical calculation up to a sufficient order in the l/m^ 
expansion. 

The current theoretical state of the art which is used in the fits is already quite elaborate. 
At leading order (which is the partonic rate) the full 0{as) and the partial 0{a1) result is 
known, while the 1/m^ and 1/m^ results are used at tree level. Hence the leading uncertainties 
are the 0{as) contributions at 1/m^ and the tree level contributions to order 1/m^. 

In the present paper we present a systematic way to perform the calculation of higher order 
terms in the l/m^ expansion. This approach is then used to perform the complete calculation of 
the 1/m^ terms for the semi-leptonic decay B X£u£, keeping the charm mass to all orders. In 
the next section we shall outline our calculational method, which amounts to a systematization 
for the tree- level calculation of terms to some order This method is then applied to the 

case of 1/m^. In section 3 we identify and discuss the parameters appearing at 1/m^, where 
we find in total five independent matrix elements, which, however, have a simple physical 
interpretation. In section 4 we discuss the impact of these additional terms on the analysis of 
semi-leptonic decays. 

2 Tree-Level Operator Product Expansion to 0{l/m^) 

In this section we outline a method which allows us - at least in principle - to calculate the 
decay rates of semi-leptonic inclusive decays at any order in the l/rUb expansion at tree level. 
The starting point for such a calculation is the hadronic tensor, which according to the optical 
theorem can be related to the discontinuity of a time-ordered product of currents across a cut. 
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Thus one starts with a correlator of two hadronic currents 

T,, = j d^a;e-"("^"-^)(5(p)|6,(a;)r^c(a;)c(0)rt6,(0)|5(p)) . (2.1) 

Here 

rM = ^7,.(l-75) (2.2) 

is the left-handed current, v — the four velocity of the decaying B meson and q the 
momentum transfer to the leptons. 

It is convenient to decompose this correlator into Lorentz scalar structure functions accord- 
ing to 

7)*!/ = -QfjivTi + v^v^T2 - ie^,uaf3v"q^T3 + q,j,quTi + (g^v^ + v^q^) Tg (2.3) 

where the scalar Tj are only functions of q^ and vq. 

Using the optical theorem we obtain for the relevant imaginary parts a similar decomposi- 
tion: 

= -g^JivWl + v^,v^W2 - ie^,ua(3v"q^W3 + q^q^Wi + {q^^v^ + v^q^) W5 , (2.4) 
and the hadronic tensor, needed for the transition amplitude, can be computed from T^i, by 

~ ImTj = Wj (left hand cut only) (2.5) 

The differential decay rate is then obtained by contracting the hadronic tensor with the leptonic 
tensor L^'^ , and one obtains 



= ^fl^ ImT^,L^^d0ps (2.6) 



where d0ps denotes the corresponding phase space element. 

Using the charged lepton energy y = 2^, the leptonic invariant mass = q'^/ml and the 
rescaled total lepton energy s = mjyVq as independent variables one obtains the triple differential 
decay rate in terms of the scalar functions Wi 



d^r _ GlmlW 



0{rnl{y^^-y^-e))e{f) 



X {w,{e, s) e + W2{q\ s) (^^^ - ^ _ + W^ie, s) e {ym, - s)^ (2.7) 

The tree-level expansion in l/m^, is most easily set up by studying the Feynman diagram 
shown in fig. 1. The double line denotes the propagator of a charm quark which is propagating 
in the background field of the soft gluons of the B meson. After rescaling the b quark momentum 
according to 

Pb — 'iTT'bV + iD (2.8) 
we write for the background field propagator 

'Sbo.^^^^ (2.9) 

where Q = rribV — q and D denotes the covariant derivative with respect to the background 
gluon field. 
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Figure 1: Tree- level Feynman diagram for the hadronic tensor in inclusive semi-leptonic B 
decays. 

The tree level OPE for semi-leptonic processes is obtained by multiplying (2.9) by the 
appropriate Dirac matrices for the left handed current (2.2). A calculation to order 
requires to expand this expression to n*'* order in the covariant derivative (iD) according to 

^^^^^ Q-rric Q - mc^^^^ Q - rric^ Q - rric^^^^ Q - rric^^^^ Q - rric^ ^ ^"^ 

We note that this keeps track of the ordering of the covariant derivatives. 

The remaining task is to evaluate the forward matrix elements of operators of the from 

kAiD^.)--iiD,Jh,p (2.11) 

where by (by) carries the spinor indices a (P). The field by is still the full QCD field, but 
redefined by a phase factor to remove the large piece of the 6-quark momentum according to 
(2.8). This field satisfies the useful relations 



b{x) 


= e-*™^^-^6^(a;) 


(2.12) 




— by i]pby 

rrib 

= ilPby + by 

2mb 


(2.13) 


P+hy 


(2.14) 


P-by 


= -^ilpby 


(2.15) 


{ivD)by 


= --^ilpilpby 
2mb 


(2.16) 



which follow from the equation of motion for the b quark field. Here P± = (1 ± ^)/2 are the 
projectors on the "large" and "small" components of a Dirac spinor. 

Note that this way of expanding will yield only matrix elements of local operators; however, 
these matrix elements contain still a nontrivial mass dependence, which will be discussed in 
the following. 

The evaluation of these matrix elements is performed most conveniently in a recursive 
fashion. The starting point is always the matrix element of highest dimension, i.e. the one with 
the maximal number of covariant derivatives. For these matrix elements one may neglect all 
contributions of order l/rUb relative to this matrix element; in other words, this matrix element 
may be treated in the static limit. 
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Thus the first step is to consider the static limit of the forward matrix element of the highest 
dimensional operator, which has the form [13] 

{B{p)\b,,a{iD^lii))....{iD^Jb,,p\B{p)) = {B,\K,a{iD^fXi))....{iD^JK,0\B,) + Oil/rrib) 

= IpaAf,^^^...^^ + sxB^^^^...^^ (2.17) 

where s\ = P^jxj^P^ is the generalization of the Pauli matrices to the case v ^ (1, 0, 0, 0) and 
|B^) is the static limit of the B meson state \B{p)). 

The tensor structures A and B have to be related to a minimal set of fundamental matrix 
elements which are defined by contracting the indices in the various possible ways. In the 
following, these matrix elements are called basic parameters for a certain order in l/m^. For 
example, at order 1/m^ these basic parameters are the Darwin term and the spin-orbit term 
Pls defined by 

2MbpI = {B{p)\K{iD^){ivD){iD'')b,\B{p)) (2.18) 
2MbpIs = {B{p)\h{iD^){ivD){iD,){-ia>''')b,\B{p)) (2.19) 

while at order 1 / ml the basic parameters are the kinetic energy parameter /t^ and the chomo- 
magnetic moment fio given by 

2MBfil = -{B{p)\K{iD)X\B{p)) (2.20) 
2MBfil = {B{p)\K{tD^){tD,){-tanbv\B{p)) (2.21) 

The hat above the quantitites means that we have defined these parameters in a covariant 
way. The definitions which have been recently used refer to the spatial components of the 
derivatives only and differ from the ones used here by higher-order terms in the 1/m expansion. 
The usual definitions of these quantities are given by 



2MBfil = -{B{p)\K{iD^)X\B{p)) (2.22) 

2MBiil = {B{p)\bSD^^){iDl){~ia^,)h,\B{p)) (2.23) 

2MbpI = {B{p)\bSD^,){ivD){iD>'^)b,\B{p)) (2.24) 

2MbpIs = {B{p)%{iD^^){ivD){iDl){-ia^,)b,\B{p)) (2.25) 

where the spatial components are 

= {g^^- - ^;V)D, (2.26) 



At 1/m^ the correction terms in the relations between the covariant and the usual definition 
do matter and are given in the next section. 

Once the tensors A and B for the matrix elements of the highest order in the 1 / mi, expansion 
have been expressed in terms of these fundamental parameters, one proceeds in a similar way 
with the matrix elements of dimension n — 1. However, now we have to take into account all 
possible Dirac structures, such that 

{B{p)\b,,^{iD,,)....{iD,„_,Xp\B{p)) = Y.^t^,,-,.-. (2-27) 

i 

where f arc the complete set of the sixteen Dirac matrices. By using the fact that the 
relations (2.12-2.16) connect different orders of the 1/mb expansion we may express the tensor 
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coefficients A^'^^ in terms of the basic parameters of the order l/m^"^ and the ones of the order 

This prescription defines a way to recursively compute the relevant matrix elements of the 
l/rrib expansion up to order at tree level, starting from the operator of highest dimension. 

Thus the leading matrix element of dimension three will then be expressed by this recursive 
method as a series in involving all the basic parameters up to this order. 

Finally, the hadronic tensor is obtained from the trace formula 

= {B{p)\KT,s^GFrlK\B{p)) = ETr{rp^^rtf« 1^(^.0) 

+ y Tr Ir^—^ Y'tt^ Y'tt^ rt f«\ ^f^i + • • • (2-28) 

where the tree level expansion of the background-field propagator (2.10) automatically yields 
the correct ordering of the covariant derivatives. Note that this saves us from computing the 
one- and even more-gluon matrix elements which would be necessary to obtain the correct 
ordering of covariant derivatives in the standard computation. 

In the follwing we explicitely perform this recursion to order and we first identify the 

basic parameters for the order 

3 Basic Parameters at 

At we have to deal with operators of dimension 7, containing four covariant derivatives. 

We find in total five basic parameters at order three of which are spin independent, while 

two are spin dependent. Written in a covariant form we define the basic parameters of order 



l/ml to be 

2Mb Si = {B{p)%iDp{ivDfiD''h„\B{p)) (3.1) 

2Mb S2 = {B{p)\KtD^{tD)HD'WB{p)) (3.2) 

2Mb S3 = {B{p)\K{{iDffK\B{p)) (3.3) 

2Mb S4 = {B{p)\KiD^{tD)HD,{-ia^'')K\B{p)) (3.4) 

2Mb S5 = {B{p)WiDpiD^iD,iDf>{-ia^'')K\B{p)) (3.5) 

which are real, since they are forward matrix elements of hermitean operators. 

In order to obtain some intuition concerning the physical meaning of these parameters, we 
may relate them to some more intuitive quantities, which are 

{E^) : Expectation value of the Chromoelectric Field squared (3.6) 

(S^) : Expectation value of the Chromomagnetic Field squared (3.7) 

{{P^Y) '■ Fourth power of the residual h quark momentum (3.8) 

((p^)(<T • B)) : Mixed Chromomagnetic Moment and res. Momentum squared (3.9) 

{{p ■ B){a ■ p)) : Mixed Chromomagnetic field and res. hehcity (3.10) 
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The basic parameters defined in (3.1-3.5) can be related to the intuitive quantities by^ 



2MbSi = -g\E^) (3.11) 

2Ms.2 = 9\{E') - {B')) + {{{pff) (3.12) 

2MbS3-(((p)^)') (3.13) 

2Mb54 = -?>g{{S ■ B){pY) + 2g{{p ■B){S- B)) (3.14) 

2MBS, = -g{{S-B){pf) (3.15) 



which also gives some information on the sign of the non-perturbative parameters. 

These parameters have to be determined independently from the lower-order parameters; 
however, for a numerical estimate it is useful to note that some of these parameters can be 
estimated in naive factorization, such as 

{B') ^ {{{S ■ B))r (3.16) 

{(pY) - {{{P'm = f^t (3.17) 

{{p'){a ■ B) r. {{pY){{S ■ B)) = (3.18) 

while the two remaining matrix elements do not have a simple interpretation in naive factoriza- 
tion. Based on the interpretation of the parameters Sj in terms of physical quantities we define 
a "guestimate" for the basic parameters Sj by 

Sl ~ f , S2 ~ -f - )Ug + , S3 ~ //^ , S4 ~ S5 /^g/^^ (3.19) 

where at least the sign of the contributions should be correctly reproduced. 

In this way all the basic matrix elements up to order 1/m^ have been identified and the 
remaining task is to express any matrix element in terms of these basic quantities according to 
(2.27). In appendix A we list all the necessary general matrix elements up to dimension seven 
in terms of the basic parameters fi,,, fiQ, Pd, Pls and S1...S5. 

Finally, we may now consider the relation between the covariant definition of the basic 
parameters and the usual ones shown in (2.22-2.25). This relation is given by 



/i^ = //^ + + ss + 455] (3.20) 

I^G-P'l- ^{Pl + Pls] - + S3 + 455] (3.21) 

pI = pI (3.22) 

pIs -pIs-^ (3.23) 



4 Results and Discussion 

The remaining task is to evaluate the scalar components of the hadronic tensor at tree level 
using (2.28) and the formulae from the appendix. We first compute the correlator (2.1) involving 

^Thc relations between the Si and the intuitive quantities are not entirely unique. For example, ((p^)^) could 
also be defined by the completely symmetrized combination of the covariant derivatives 
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the time-ordered product. From the scalar components shown in (2.3) we need only Ti, T2 and 
T3 due to current conservation of the leptonic current. The resulting expressions up to order 
are tedious and are given in appendix B 
Taking the imaginary part to obtain the components of W^i, according to (2.5) we use the 
relation 

TT n! 

from which we can obtain the triple differential rate up to at tree level using (2.7). The 

resulting expressions of the double differential rate, the single differential rate and the total 
rate are quite lengthy and given completely in appendix C. 

The relevant quantities for the experimental analysis and the determination of Vet, ^-re the 
moments of the lepton energy spectrum and the moments of the hadronic invariant mass spec- 
trum. In order to get a quantitative idea of the effect of the terms we consider the 
to the moments, normalized to the partonic rate at tree level. Thus we define 

6i'HM',)^-jm.M', J mj^^ (4.2) 

Ecut 

,i'\Et)^-jm. J iE,Et^^^ (4.3) 



cut 



To = ^^^i^i^ - 8P - 121og(p)p^ + - p') (4.4) 

where F*^^) is the contribution of order l/m\. 

We first study the dependence of the 0{l/ml) contributions to the moments on the different 
non-perturbative parameters Sj. We write the moments as 



5(^)(Ml) = ^m^i^(")(^,,0-^ (4.5) 
5^'\E^) = j^mlgr\E^^,)^^ (4.6) 



i^ 


1 


2 


3 


4 


1 


- 11.570 


- 6.314 


- 3.418 


- 1.845 


2 


2.073 


1.076 


0.556 


0.288 


3 


- 5.969 


- 2.801 


- 1.320 


- 0.620 


4 


- 0.102 


- 0.126 


- 0.105 


- 0.074 


5 


- 3.377 


- 1.042 


- 0.174 


0.089 



Table 1: Values of the coefficients gi with a weak dependence on £^cut- The values quoted are 
for Ecut = 0.8 GeV. 

The contributions of order 1/m^ are strongly concentrated in the endpoint ~ (m^ — 
ml) I (2mfe) for the lepton energy. Hence the dependence on the cut-off energy of the coefficients 
is small and can be neglected, except for the functions jf^ ^ corresponding to the first hadronic 
mass moment. 
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i\" 


2 


3 


4 


1 


6.214 


- 6.633 


- 1.322 


2 


- 1.343 


1.026 


0.203 


3 


2.472 


1.358 


- 0.033 


4 


- 0.059 


0.315 


0.133 


5 


- 0.377 


- 0.129 


0.019 



Table 2: Values of the coefficients f-^^ with a weak dependence on -Ecut- The values quoted are 
for Ecut = 0.8 GeV. 
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Figure 2: The dependence of fj;^^ on the energy cut in the lepton energy. 
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Figure 3: The dependence of f-^^ on the energy cut in the lepton energy, leaving out the much 
larger function fl^\ 
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In tables 1 and 2 we tabulate the values of the coefficients in (4.5) and (4.6) for those 
coefficients which are practically independent of -Ecut- 

As pointed out above, only the coefficients f-^^ of the first hadronic moment depend on E^nt 
in a substantial way. Figs. 2 and 3 show the dependence of the five functions on the energy cut. 

Finally we shall also study the overall effect of the 0{l/ml) contributions. In order to do this 
we have to insert numerical values for the nonperturbative parameters. We use the relations 
(3.19) and the values obtained from the fit in [10] for the lower order basic parameters. In 
table 3 wc list the values from [10] for reference; in table 4 we list the values of the basic 
parameters at 1/m^ obtained from our guestimate. 



Parameter 


(GeV)2 




Pis (GeV)-^ 


pI {GcVf 


Value 


0.401 


0.297 


-0.183 


0.174 



Table 3: Values taken from Buchmiiller and Flacher. We do not show the uncertainties, which 
are at the level of 15%. 



Parameter 


si (GeV)^ 


S2 {GeVy 


S3 (GeV)^ 


S4 (GeV)^ 


S5 (GeV)4 


Value 


-0.076 


0.148 


0.161 


-0.119 


-0.119 



Table 4: Our guess for the basic Parameters Sj which is used in the numerical analysis. 



Inserting for the masses the values = 4.59 GeV and = 1.142 GeV from [10] we obtain 
for the contributions to the moments for a lepton energy cut of 0.8 GeV the values shown 
in table 5. As pointed out above, the I /ml contributions to the lepton energy moments are 
practically independent of Ecut, while the dependence on Ecut of the hadronic moments is given 
by the functions f-^^ . 



n 


1 


2 


3 


4 




-0.1835 GeV 


-0.0104 GeV^ 


0.185 GeV=^ 


0.1064 GeV^ 


5W{E2) 


0.0066 GeV 


0.0154 GeV^ 


0.0351 GeV=^ 


0.0803 GeV^ 



Table 5: Numerical values for the contribution of order 1/m^ using the parameters from table 4. 



The contribution to the total rate is very small; using our estimates for the parameters Sj 
we get 




resulting in a completely negligible shift of the central value of Vcb from the terms of the order 
l/ml 



5 Conclusions 

In this paper we presented the complete result for the 0{l/ml) contributions to the semileptonic 
rate at tree level. Although wc have in total five new, non-perturbative parameters, we can use 
the known matrix elements to order 0{l/ml) to estimate the size of the contributions. 

It turns out that the size of the 1/m^ terms is "normal", i.e. we do not see any abnormally 
large coefficients in the l/m^ expansion, at least in this tree level calculation. 
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The formulae collected in the appendix allow us to include now the Oil/mf) contributions 
into the moment analyses of semi-leptonic decays. In particular, the fourth central moment 
in M\ is dominated by the terms, so these terms are required for a sensible analysis of 

higher-order moments. 

The method suggested here allows in principle the calculation of arbitrarily high orders in 
the l/m;, expansion at tree level. However, the limitation is the number of basic parameters, 
which is presumably growing factorially, leading to a factorially growing effort to express the 
general matrix elements through the basic parameters. 

On the other hand, assuming that the ^/fii^^ terms also behave "normally", the uncertainty 
from these contributions in the determination of Vcb will really be extremely small. The un- 
certainty currently assigned to the use of the HQE for the detemination of Vch is about 1%; a 
conservative estimate is to use the size of the last term calculated as the uncertainty means 
that the uncertainty in the Vcb determination due to HQE may be reduced to < 1%. 
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A General matrix elements in terms of basic parameters 

In this appendix we list the general matrix elements up to dimension seven in terms of the 
parameters /i^ (2.20), fic (2.21), pd (2.18), pls (2.19), and S1...S5 (3.1-3.5). The results are 
written as Dirac matrices which have to be inserted into the trace formula (2.28). MATHEMATICA 
notebooks with the corresponding expressions can be obtained from the authors. 

A.l Dimension seven 



{B{p)\ UtDP){iD''){iD^){iD')K\B{p) ) 

Mb, 
30 



+ i-ia^^) {g"" - v^v") ^ (4s5 - S4) 



Mb 
+— P_ 
60 



+ {-laY" {g^'' - i;V) (2^4 - 3^5) |p+ 
{v^v^ - g^^) {gP" - v^v") (-2si + 3^2 - Tsg) 



+ 2 {g^" - v^") {g 



') (Sl + S2 + S3) 



+ {gP^ - v^yP) {g^^ (-8si + 7s2 - Ss^) + v^" (28si - 7s2 + Ssa)) | 
+ 0{l/mb) 



A. 2 Dimension six 



(A.l) 



{B{p)\ K{tDP){tD-){tD')b,\B{p)) 
Mb, 



3 

Mb 



6 

Mb 



6mb 

Mb 



P4 
P. 



12mb 

Mb 



24mfe 
Mb 



Mb , 



24m6 



+ 



V v 



{{-iaP'')v,v\'' - {-ia^^)v,vPv'') (3si -82 + 33-3^ + S5) | 
{{-iaP'')v,v^ - (ia^^HvP) r (-2si + ^2 - S3 + Sss) | 

((_i(jP-)^A _ (_i^Aa)^p^ ^_2g^ + S2-SS + 2S5) |p+ 



12 



Mb 
24m5 

Mb 
60mf, 

Mb 
UOrrib 



+ 



Mb 
Urrib 

Mb 



Mb 
Mb, 



Mb 
Mb 

Mb 
'uOrUb 

Mb 
'uOrrib 

Mb 

12mb 

Mb 



((■ 



la' )v — [—ta jv 



(7^ - V^f) (Si + S2 + S3 + S4 + S5) 

(^P^ - v^') ((-6P+ - 2)si + (4P+ + 3)s2 - (7 - 4P+)S3 

-10S5+^(8P+ + 6) (S4 + S5))^"'^| 



(7V - v^) (^(8 - 6P+)si + (4P_ 
+ (4P+ + 3) (^7^-^;^) sa)^'"^! 



3)S3 + IOS5 - -(14 - 8P+) (S4 + S5) 



+ A"'^) (S3 + 55) 



v^'t;''7'^ (si + S2 + S3 + S4 + S5) 



-2si + S2 - S3 + 2S5) 



V (33si - 7s2 + 3s3 - IOS5 - 2 (s4 + S5)) 

^;P^;^^A ^ 3^^ _ ^ ;^Q^^ - 2 (^4 + S5)) | 

v^fv^f (-12si + 3s2 + 13s3 + 3s4 + 13s5)) I 
+ 0(l/m^) 



(A.2) 
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Dimension five 



{B{P) I 
Mb 

Mb 
48ml 

Mb 
ASml 
Mb 
6m6 
Mb^ 

Mb 

Mb 
'l2ml 

Mb 



K{iDP){iD'')K\B{p))^ 

( nr.P'^\ i 9r,2 , -2^2 - g4 + g5 2 (pf, + p3g) 
^"^^ + 

'-io^") (s2 - S3 + S4 - Ssg) 

/3S2 + 3S3 + I6S5 8(p|5 + p^g) 



I- 



% rrib 

{v^Y + ^"^7") (4si - S2 - S3 - 455) 



(Sl - S2 - S4) 



v^t;'^^ (-12si + 3s2 + S3 + 2s4 + 4s5) 
-rO{l/ml) 



Dimension four 

{B{p)UiD^)K\B{p)) = - ^P^[vM^l - 1^1)] + ^{(7" - 

Dimension 3 

{B{p)\KK\B{p)) = P+Mb + ^{ftl - f^l) + 0{l/ml) 
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B Scalar Components of the Time-Ordered Product 

To compute the differential rate in tlie limit of vanishing lepton masses only the functions Ti, 
T2 and T3 are needed. In the following subsections we give the complete expressions for these 
scalar functions up to order We use the notations 

Aq = ml — 2mt,q ■ v — ml + (f (B-l) 

Furthermore we use dimensionless variables according to 

2 ^ ^0 A mo^ 

— ^ q , mb- V ■ q ^ mbV ■ Q , — ^ ^ Ai (B.2j 

m^ mi mi 

MATHEMATICA notebooks with the corresponding expressions can be obtained from the authors. 



B.l Scalar Function Ti 

v-QMb 



Mb.2 f _ 1 V- q{5v ■q-3)-2q'^ _4vQ{{v- qf - q^) 



m\ 



6A1 3Af 3Af 



Mb .2 ( ^ -llg' + 10(t;-g)' + 7p+l 
+ Smg^^ Ui ^ 2A2 

Mb .3 / 2 4(vg-g2) + 3p + 3 A{{vqf -q^)vQ 

"I tPd 77a ^ 



ml 



+ 



3Ai SA? 3Af 



3A^ 



Mb .3 [2 -Aq^ + A{vqf + Zp-l Av Q {{v qf - q^ 



+ —A Pis ^ + 



^4 f^Ls y^^^ 3^2 3^3 

Mb (2vQ 2 (-30(^; • + 38(^; • + 15 - 4) ^; • ? + + 30) 



8 (g^ - • qf) (-6g^ + 2^ • q{?>v • g + 5) + 5(p - 2)) 

15Af 

16 {{v ■ qf - q^) (-6g2 + 6(^; • g)^ + 5p) v • Q 
15Af 

Mb /l-3w-g (35^ • g - 23)g2 - 2v • g(?7 ■ 9(25^ • g - 39) + 40) + 40 
^7^^' V 6A2 + 30Af 

2^2 



2 (g2 - 6(v • qf + 5p + 5) (g^ - {v ■ qf) ^8v-Q{q^ -{v- qf) 



15At 5Af 
Mb /-y-g + S 3(5v • g - l)g2 - 2^ • g(v • g(15w • g + 11) - 40) - 40 



^ V 6A2 30Af 

4 (3g2 + (10 - 13v ■ q)v ■ q) (g^ - {v ■ q)^) ^8v-Q{q^ - {v- qff^ 



15Af 5Af 
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Mb ( _ 7vq + l 2{vq + l) {q^ - {v ■ qf) 



ml^\ 12Af ' 3Af 

2{q^-{vqf) {q^-2{vqf + p+l) \ 

Mb f3v-q + 5 2 {q'^{3v ■ q + 1) - 2v ■ q{2 - 3v Qv ■ q)) 



ml 



4AI 3AI 
, 4 (g^ - (t; ■ qf) (-1 V + 13(t; ■ q)' + 16p + 6) ' , 
+ 15A^ ' ^^-^^ 



B.2 The scalar function T2 

2Mb 



2 



Mb. 2 ( 5 14wg 8(g2-(y.g)2) 



m^^" V 3Ai 3A? ^ 3A? 
Mb .2 f 5 A-lOvq 



3mr'' V ^1 ^? 

Mb ( 4 4-4?;-g 16f • • g 16t' • Q (y- - (r • fy)^) 



3Ai A? 3A? 3At 



Mb f 4 2-3v-g 4 (g^ _ 3^ . . 



ml 



3Ai 12A2 3A? 



Mb f A{3v-q + b) 4 (-39g2 + 42(v • g)^ + 40^; • g + 35p - 5) 



ml ' \ 3A2 15Af 

32 ((6v • g + 5) ((w ■ g)2 - q^) + 5v • gp) 
15Af 

32 (g2 - (t; . qf) (-6g2 + 6(v • qf + 5p) 



+ 



15Af 



Mb f5v-q + 3 -39q'^ + 62{v ■ qf + 5p + 15 
?r^2 — 77T9 ^ 



mf V 3A? 15A? 



2^2 



^ 8(6v • g + 5) ((v • g)2 - g2) ^ 15 (g^ _ . 



15Af 5Af 
Mb /3i; • g - 1 + 82(t; • qf - 25p + 5 



^ mf V 3Af ^ 15Af 



, 128T;-g((^-g)^-g^) , 16 (g^ - • g) 



2^2 



15Af 5Af 
Mb f Av-q + 1 4 (g^ + (5 - 8v • g)^ • g) 4^ • Q (g^ + 2(1; • g)^) 



V 6 A? 15Af 5A^ 

Mb /24f ■ g + 13 -94g2 + 172(t; ■ qf - 2Qp + 20 



^5 TTTo. + 



ml ^ \ 6Af 15Af 



, 4 ((17 - i7v ■ q)q^ + 2{v ■ qf {I3v q + 2)) \ 

+ i5Al ' ^^-^^ 
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B.3 The scalar function T3 



J-2, — 



Mb ^2 f G — 5v ■ q 



Mb .3 f6-4vq Av ■ Qv ■ q _ 8v ■ Q (q^ - {v ■ qy)\ 
^7^^^ V 3A2 + ~3Af 3A| ; 

+ V 3A? ^Af~ ; 

Mb (2 (g2 - 6(t; ■ g)^ + • g - 6) 16(3^; • g + 5) (g^ - • qf) - AQv ■ qp 
^7^^' V 3Af ^ 15Af 

16 (g2 - (^; . qf) (-Qq^ + 6(v • qf + 5p)~ 
^ 15Af 
Mb f 5 -9q^ + 10{v ■ qf + Ap + 2 
^7^^' V3Af^ 6Af 

4(3^; ■q + 5){{v qf - q^) ^{q^ -{vqff \ 
lbA\ ^ 5Af J 



Mb ( 1 5g^ + 6((^-g)^-p) 52t; ■ g ((t; ■ g)^ - g^) 8 (g^ - (t; ■ g)^)^ 

^ mf^ \ Ay 6A? ^ 15Af ^ 5Af 

Mb /I 2 

H frS4 



(p- (f -g)^) 4^-g((^;-g)^ -g^) ' 

QA3 CA4 



m^' V4Af 3Af ^ 5Af 

Mb / 35 



+ — S5 



, , llg^ + 12(r-g)^ + 9p-9 , l(13r • g - 25) ((r • g)^ - g^) ^ 

l,12Af + 3Af + 15A| i ^^-^^ 
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C Decay Rates 

Finally we list in this appendix the full expressions for the differential and the total decay rate 
at tree level up to order I /ml. MATHEMATICA notebooks with the corresponding expressions can 
be obtained from the authors. 

C.l Double Differential rate d^T / {dEo dM]^) 
We use for this double differential rate the variables 

Eo^rrib-v-q, ^ {mi,v - qf (C.l) 

£ =^l^-l^ ^^^^A (i6 {aeI - 3 [mi + 1) ^0 + 2Mi) - p) 

+ 3^A' ((404'-2io'(lOMi + 5p+ll) -3p 

+ £^0 (ep - lOM^ ) + 5Ml + 7{p + l)Mi^ S'{M], - p) 

+ ^16£^^ - 4 {2Ml + P + 3) ^0 - 8P^0 + [m^ - 5p) 
+ (-Ml + (13p + 3)Mi + 9p) 5"(Ml - p)^ 

+ ^/iG (- ( + 40Eo3-2^o'(lOMi + 5p+ll) 

+ io (6p - lOMi) + 5M| + 7(p + 1)M2 - 3p J <5'(M|- - p) 



{%El - 3 (Mi + 1) io - 2Mi) [mI - p) j 5"(Mi - p) 
- 2 (^16^^ - (4Mi + 2p + 5) + (ep - 4Mi) 
+ 2Ml + (p + 4)M|. - 3p^ (5'(M|. - p) 

-4(E2_M2j (442„ (2m2 ^^^;^^^^_^2 5"(Mi-p) 

- ^ (^32^^ - 8 (2M|. + p + 3) + 8 (p - 3Mi) 

+ 2 (SM^ + 2(5p + 6)Mi - 3p2) ^0 + 2 (sM^ - 22pM| + Qp^) 



- 3 (mi + 3) (Mi -p) Y'\Ml-p) 

+ 3^pi5 ( - 2 + 164' - 24' (4Ml + 2p + 5) + 4 (6p - 4Mi) 
+ 2M| + (p + 4)Mi - 3p ) (5'(M2 - p) 
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-4(io'-^l) (4io'- (2Ml + p+l)io-Mi + p) 5"{Ml-p) 
+ ^ {sEl (mI - p) - 2£:o' (^x + 2Mi - p(p + 2)) 

- 2Eo (^iV4 + 2pMl - 3p2) - Mi + (6p + l)Mi 

+ (2 - 5p)pMi - 3p2 V^'H^i - P)^ 



4 

+ 



2 



3ml V V 

+ (qM^ + (3p + 8)Mi - 3p) io - 2pMi^ 5"(M|: - p) 

- ^ - (^2AEl - 6 (2Mi + p + 1) 4' 
(Sp - 9Mi) 4 + Mi (2Mi + p + 1) ^ (M^ - p) 

^ (^1924' - 48 (2Mjr + p + 3) 4' + (48p - 176M|.) 

8 (^9Mi + (17p + 21)Mi - 5p2) ^0 + 4 (l9Mi - 64pMi + 25p2^ E^ 

- (^21M|. + (69 - 2p)Mi + 5(p - 18)pM|. + 45p2) ^0 

- 2Mi (^M^ - lApMl + V) - P)) 

+ 3^52 (^4 (^40£'o" - 2 (lOMi + 5p + I3) El - 4 (m^ + 2p) 4' 



+ 
+ 



+ 15 



+ (2M^ + (19p + 8)Mi + 9p) Eq - lOpMl j 6"{Ml - p) 

^ (^48io' - 12^0' (2MI + p + 1) - 6^0' (llMi + 5p) 
2£^o' (l4Mi + (17p + 7)Mi + 5p(p + 1)) + 6£^o (sM^ + 5p2) 
4Mi - (7p + 2)iV4 + 5(1 - 5p)pM| - ISp^^ 5^'^\mI - p) 



+ ^ (^192^^ - 48 + p + 3) 4' - 16 (llMi + 7p) E^ 

+ 24 (^3M| + (9p + 7)Ml + 5p) + 4 (^19Mi - 34pM|. + ISp^) El 

- 3 {lM% + (26p + 23)Mi- + 5p(3p + 2)M|. + ISp^) Eq 

+ 3^53 (4 (^AEl - 2 (6Ml + 3p + 7) El + (8p - 12Mi,) El 
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+ 



2 

+ 



3ml 



+ (gM^ + (3p + 8)M2 - 3p) Eo - 2pM|^ 6"{M]^ - p) 

+ ^ (^208i^ - (^26M^ + 13p + Ss) - 2El (^23M|. + 65p) 
+ 2El (lAM"^ + (67p + 27)M| + 15p(p + 3)) - 2Eq (21M^ - 5QpM]^ + 45^^) 
+ 16Afj. + (18 - 67 p)M^ + 15(p - 5)pM| + 45p^^(5(^)(M|. - p) 

^ (^192^^ - 48 (2M^ + p + 3^El-16 (m^ + 17p) E^ 

16 (^2M^ + (lip + 3)M^ + 5p(p + 3)) - 4 (^21M| - 26pM^ + Sp^) 
+ (^19M|- + (51 - 38p)M^ - 25p(5p + 6)M|. - 45p2) £"0 

- 2Ml (^M^ + 6pM| - 55p2) ^ ^^^^ (M^ - p) 
S4 (^32E^ - 8 {2MI + P + 3) £3 + 8 (^M^ - 3p^ El 

+ (-5M1 + (29p + 3)Ml + 21p) Eq + 2Mi (^M^ - 9p] ^ 5"(Mi - p) 
+ ^ (^244^ - 6^^ {2MI + p + 1) + 2El {mI - 25p) 
+ El [-6Mx + (47p - 13)Ml - 5(p - 5)p) + 2^ (2M^ - hpM^ + 1 V) 
+ 3M% + (4 - llp)Mi + 5(1 - 2p)pM|: - lbp^^5^^\Ml - p) 

^ 'eI - Ml) [mI - p) [-&EI + 5 (ifi - p) + Ml + 5p) 5^'\Ml - p) 



+ 



Y — — ^ J Y — / \ \ — / — — / y I 

(^192^0^ - 16 (6Mi + 3p + 7) El - 8 (l3Mi - 9p) El 
+ (^51Ml + (21p + 67) Mi - 27p) ^0 + 2Ml (m^ - 9p) ^ 5"{Ml - p) 
+ ^ ^104£^^ - 2El {2QMI + 13p + 33) - SEl (iQMl - 5p) 
+ El {lAMl + (87 - 23p)Ml + 15p(3p - 1)) + Eq [-QMx + 80pM| - 90p2) 

- 7M% - (lip + 6)M^ - 45pMi + 45p^^ 5^^\mI - p) 

- ^ (^108 (mI - p) Et - 30 {2MI - 3pMl + M|. + (p - l)p) El 

+ (-mMl + 16pMl + 70p2) El + 15 (^3Mi + (1 - 4p)M| + p^M^ - p^ j 

+ 8Ml (Mi + 4pMl - 5p2) ) (M^ - p)) | (C.2) 
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C.2 Differential rate dV/dy 

We use in this Differental rate the dimensionless variable 

y = 2^ (C.3) 
rub 

Furthermore x is a short hand notation for x = y^(1 ~ V ~ p)- means therefore the nth 

derivative of the delta function with respect to its argument x. To rewrite the argument and 
derivative of the delta function linear in y, we use 



where 



1 dy _ Id Id 

y{y) dy) g'{y) dy " ' g'{y) dy 

" V ' 

n times 



dy - 1927r3 '^^^1 \ {y - 1)3 [y - 1)2 2y + 3J 

2y^ .2 /2(y2_,5y + 10)p3 (15 - 6y)p2 



2y2^2 /^ 5(t/^-% + 6)p3 9(^^)p2 g(2y-3)p 

^^^^ I ^^Tj^ {y-ir ^ 

2 .3 (5y^ - 30y3 + 75y2 - 84y + 54) p3 

^ yM-V + 25|/^-50|/ + 42)p^ ^ 2y2 (2y2 - 3?/ - 3) p 



^ 2y^2y' - 5y + 9) (p + l)(p + 2)(p - 1)^ 
y - 1 p2 



p4(V + 3y + p-2) p^y , (y-2)p2 p 



3m3 



V (y-1)^ y-ir ^ ^ 

.3 / (y3-5y^ + 10y + 6)p3 3(^2,4^^6)^2 6(2y - 3)p , ^ \ 



1 / _ %^ (1 V - 11 V + 357y3 - 595y2 + 520y - 240) 



15mt 'V (1/ - 1)' 

2y2 (43y^ - 258y3 + 645^^ - 760y + 420) 



+ 



{y-ir 

20y2 (4y3 _ 15^2 ^ 20y - 12) p AOy^ (2^^ - 4y + 3) 



+ (y-l)'i + (y-l)2 

2 (3p^ + 10p3 - I4p2 - 34p - 45) {p - if ^ ^ _ 

p3 
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^ _ 3(|/-19)p5 ^ (-2y2 - 28y + 75) 2{by + 8)p3y 



{y-l) 

2{lly + 2)p^y 3 (y^ + 4y + 5) p _ \ 



6 



+ + 



(y - 1)4 (y - 1)2 

3p6 (5-2y)p5 2p4^ 



(^_l)io (^_l)9 (^_l)7 

, (3y- 5)p2 3p , 



(y-l)5 (y-l)3 y-l 
1 /^6y2 (27yS - ISQy^ + hQly^ - 945^^ + 820y - 340) p^ 

+ — I'S2 



60m^ V (l/ - 1)^ 

6y2 (21y4 - 126^^ + 315^^ - 400y + 220) p^ 40y'^ {by^ - 16y + 14) p 

, 40(y-2)y^ 2 (3p^ + 23p^ - lip + 45) (p - 1)\ , ^ 

+ -(7rTj^ % + 

/3(y - 19) p^ (12y2 + I28y - 185) p^ _ 2 (15y2 + 32y - 105) 

+ V {y-iy ^ (^^1? 

2(6y^-8y-45)p^ (3y^ + 12y + 5) p \ 
^ , 3p6 , 3(4y-5)p5 , 6(5-3y)p4 



(y-l)io (y-l)9 {y-iy 
^ 6(2y - 5)p^ _ 3(y - 5)p^ _ ^,^„^^^ 



(y-l)5 (y-l)3 y 

2y2 (-29y^ + 203^^ - 609^^ + 1015^^ - I240y + 480) p^ 



2y2 (7y4 _ 42y3 ^ io5y2 - 400y + 240) p^ 
_ 40y^ (3y=^ - lly^ + 14y - 3) p _ 40y^ (4y^ - 14y + 9) 

(7^1)^ 

^2(17p^ + 34p + 45)(p-l)^ 

p3 

^ (37y + 17)p5 ^ (52y^-32y-65)p4 ^ 2 (5y^ - 12y + 45) p^ 



{y-lf {y-lf {y-lf 

2(14y^-32y + 25)p^ (3y^ - 28y + 5) p „\ ^ 

/ 3p6 3(4y-5)p5 6(5-3y)p4 
V 

6(2y - 5)p^ _ 3(y - 5)p^ _ 

+ (y-l)5 (y-l)3 y-lj^^^^^ 
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+ 



10ml 



S4 



2y^ - 6y^ + 15y^ - 20y - 10) _ y^ {7y^ - 35y^ + 55y + 45) 
2y^2y^-5y^ + 15)p 



(y - 1)^ 



+ 



(y - 1)^ 

2p5 



+ 5y^ + 



2{p-2){p+l){p-iy 



■6{y + p-l) 



6{y-2)p^ 2{3y^-13y + 12)p' 



(y - If (y - ly 
2(y2_7y+10)p2^2(y-3)p 



(y - 1)^ 



+ 



S5 



4|/2 (19|/^ - 114|/3 + 285y2 - SSOy + 60) p^ 
^ y2 {-lly^ + 55y2 - 65y + 285) ^ 2|/2 {7Ay^ - 215y^ + 120y + 165) p 



+ 



{y - 1)^ 

15(19 - lly)y2 2(13p2-3p + 74)(p-l)4 



y-1 



+ 



(y-i)' 

6{y + p-l) 



+ - 



26p5 ^ 2(19?/-8)p^ ^ 2(y2-3i^/ + 54)p3 



(y-iy 



2 (7y2 - 39y + 80) p2 2(31 - 7y)p 



{y - 1)^ 

C.3 Total decay rate F 



+ 



(y - 1)^ 



(y-i)^ 
y'5'(a;) 



(C.4) 



1927r3 
1 



2m 



+ 



+ 



+ 



+ 



2m 



3ml 



9ml 



9ml 



-2 



S2 



4 



+ 



9ml 
1 

3ml 



S5 



p^ + 8p3-121og(p)p2-8p+l 
p^ + 8p=^-121og(p)p'-8p+l 



- 5p' + 24p-^ - 12(log(p) + 2)p' + 8p - 3 



- 5p^ + 16p=^ - 12p2 - 16p + 12 log(p) + 17 



9p^ - 20p^ + 9p2 + 6 log(p) + 2 



- 27p^ + 76p^ - 72p2 + 36p - 12 log(p) - 13 



3p' - 7p2 + 9p2 - 21p + 4(3 log(p) + 4) 



- 5p^ + 16p^ - 12p2 - 16p + 12 log(p) + 17 



(C.5) 
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